Symmetries of differential-difference dynamical systems in a 

two-dimensional lattice 

Isabelle Ste-Marie* and Sebastien Tremblay"'' 
Departement de mathematiques et d'informatique, 
Universite du Quebec, Trois- Rivieres, Quebec, 
G9A 5H7, Canada 

June 29, 2009 



Abstract 

Classification of differential-diiTerence equation of the form Unm ~ -Fnm(i, {^'pg}|(p,q)gr) 
are considered according to tlieir Lie point symmetry groups. The set F represents the point 
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dimensional for nonsolvable symmetry algebras. 
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1 Introduction 



The purpose of this article is to perform a symmetry analysis of the following class of differential- 
difference equations 

^nm = Unm — Fnm{t, {Up?} I (p,(j)er) = 0) (1) 

where the overdots denote time derivative. The set F being the point (n, m) and its six nearest 
neighbors in a two-dimensional triangular lattice (see Figure 1) given by 



r=<.{n,m), {n+l,m), {n,m+l), (n— l,m-|-l), (n— l,m), (n, m— 1), {n+l,m—l 




The dependent variable Unrn{t) can be interpreted as atomic normal displacement from their 
equilibrium position to the 2-dimcnsional triangular lattice, where (n, m.) locate the vertex on skew 
coordinate system. The function Fnm, called the interaction, is an a priori unspecified smooth 
functions. Our aim is to classify such system according to the Lie point symmetries that it allows, 
that is, to classify the functions Fnm- 



(„_!,„ + !) {„,m + l) 



• • • 

(n-l,m) (»,™) (n+l,n.) 



• • • • 

{n,m-l) {„+l,„_l) 




Figure 1. The atoms (n, m) and its six neighbors in a two-dimensional triangular lattice. 

The assumption for this model are the following: 

(i) The interaction Fnm involves only nearest neighbors in the triangular lattice, i.e. atom (n, m) 
interacts only with the atoms of the set F, see Figure 1. 

(ii) In the bulk of the article the interaction Fnm is assumed to be nonlinear and coupled, i.e. 

8'^F f)F 

" ""^ 7^ for some (p, q), {p', q') G F and ^ ^ for all {p, q) G F. 

(iii) The interaction Fnm is isotropic in the three directions of the lattice. 

(iv) We suppose that interaction Fnm depends continuously on the discrete variables n and m. 
For instance, interactions depending on terms of the type (—1)'"" are excluded. 



(v) Only "maximal" symmetry algebras for a given interaction F„m are listed. In other words, if 
a given interaction Fnm allows symmetry algebras £i, . . . ,£n with dim£i < dim £2 < • • • < 
dim>Cjv, then wc will only list the case Cn- 

Our motivation is the same as for classifying differential equations according to their symmetries, 
see [UlllinJU- When (P) allows a nontrivial symmetry group, then it is usually possible to obtain 
exact analytical solutions satisfying certain symmetry requirements. Model of this type have many 
applications in solid-state physics, see for instance the work of Biittner et al [SJ [51 [3 [H] . 

The formalism used in this article was called "intrinsic method" [H [10] . It has already been 
applied by Winternitz et al to the following one-dimensional cases: monoatomic molecular chains 
[TT] . diatomic molecular chains [T^ and to a model with two types of atoms distributed along a 
double chain [13]. In this paper we consider a two-dimensional case. In this formalism the Lie 
algebra of the symmetry group, often called "symmetry algebra" , is realized by the following vector 
fields 

X = T{t,Unm)dt + <f)nm{t,Un,n)du^^- (2) 

The algorithm for finding the functions r and (j)nm is to find the second prolongation pr^^^X of X 
and to impose that it should annihilate ([T]) on their solution set, i.e. 



pr 



= 0- 

A„„=0 



Our prime objective will be to find and classify all interactions Fnm for which ([T]) allows at 
least a one dimensional symmetry algebra. Thereafter, we will specify interactions further and find 
all those that allow a higher dimensional symmetry algebra. 

Note that there are other formalisms to find symmetry algebra of differential-difference equa- 
tions. Consequently the use of the intrinsic method - applied for fixed lattice and used here - is 
not obligatory. Another method, called the "differential equation approach to differential-difference 
equation" J^, can also be considered for such kind of lattices. Other methods exist in which the 
group transformations can also act on the lattice [Ml [151 IHl [12] • A complete review of these 
different methods have been considered in IS]. 

In Section 2 we present the determining equations for the symmetries and illustrate the allowed 
transformation, i.e. the classification group that will be considered to obtain class of equations. 
In Section 3 all the calculations to find the one-dimensional symmetry algebras are explicitly 
given. Some notations are also introduced for the rest of the paper. Section 4 is devoted to 
interactions with abclian symmetry algebras. They are denoted by Ajk, where the first index shows 
the dimension of the algebra and the second index enumerates the nonequivalent classes. Their 
dimensions satisfy 1 < dim£ < 12 with dim£ ^ 9,11. Nonsolvable symmetry algebras, denoted 
by NSjk, are treated in Section 5. They all contain sl(2,]R) as subalgebra and 3 < dim£ < 13 
with dim£ 7^ 10, 12. The results are summed up in Section 6 where we also outline future work 
to be done. 

2 Formulation of the problem 

The second prolongation of the vector field (Hj) is given by [9] [TOl [11] 
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(p,9)er 

The coefficient is a function depending on n, m, Unm, "nm and u„m given by 

'/'nm = Df<j)nm ~ {DIt) Unm " '^{Dit) Unm, (5) 

where D( represents the total time derivative. 

From ([3]), ^ and ([5]) one can obtain the determining equations for tlie symmetries. We 
eliminate the ilnm terms using ([T]) and then request that the coefficients of u^^, Unm and 
should vanish independently. From the determining equations of the coefficients k — 1,2,3, 
the vector field ^ must have the form 



X = T{t)dt 



(6) 



where a„jn — 0. The remaining determining equation involves explicitly the interaction Fnm and 
is given by 

■f •• / 3 \ 



E 

(p,«)er 



(7) 

^Ujjn ^ nm ■ 



Our aim is to solve ([7|) with respect to both the form of the nonlinear equation and the symmetry 
field. For every nonlinear interaction Fnm we wish to find the corresponding maximal symmetry 
group. Since for any symmetry group there will be a whole class of nonlinear differential-difference 
equation related to each other by point transformations, we will look for the simplest element of 
a given class. Hence, we shall classify ([TJ into equivalence classes under the action of a group of 
"allowed transformations" . We restrict the allowed transformations to be fiber preserving, i.e. to 
have the form 

i=i{t), Unrnit) = ^nm(t,Unm{i)j, {n, fh) = {u, m) , 

where flnm and t are some locally smooth and monotonous (invertible) functions. Substituting 
these transformations into ([T]) and requiring that the form of the equation is preserved, we find 

^nm(^) — t ^ Pnm ^nm(t^ ~^ Qnm{t) 7 

where Pnm 7^ 0, Pnm — and < 7^ 0. The transformed system is then given by 

y^nmii) — Fnmiij {^P<?} I (p,g)(=r) ; 

where 

Fnm i ^nm {^nm ^Tim(^)^ ^^^'^^ ^ ^ ^nm{t^ • 

The vector field given by ([6]) is transformed into 
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X = 



-{t)idl + I + flnm + \i ^ i Unni 

In other words, under the allowed transformations the vector field © characterized by the triplet 
{r(t), Gnm, A„m(i)} is transformed into the triplet {T(i), a„m, ^nm{i)} where 



■-{i) = T{t{i)) i, 



3 One-dimensional symmetry algebras of the systems 

Let us now introduce some notations that will be useful in what follows. First we will write 
functions g of the type 9{{£,pq}\{p,q)£s) j where S" is a subset of and : I? ^ M., more simply 

as gi^pq) with (p, q) e S. Second, if {fpl\ fpq\ . . . , fp^^} is a set of A'' functions fpq : ^ M 
depending on the discrete variables p, q, then we define the determinant function T) : 1?^ ^ M by 



(2) 

P2<32 ' 



r(N) 



PjvJnJ 



Jpiqi 



Jp2q2 



A2) A2) 

Jpiqi JP2q2 



AN) AN) 
Jpiqi 



^(1) 

JpNqN 
f(2) 

jPNqN 



AN) 
JpNqN 



For example, we have I?[/„„i, g„+im, l«m+i 



fnni fn+lm fnm+1 
Qnra Qn+lm gnm+1 



where the function 1 



pq 



^ K is the constant function (p, q) i-^ 1. (Here, the indices of the constant function are written 
only for clarity of the determinant function 'D[fnm, gn+im, Inm+i]-) 

Theorem 3.1. Equation ^ allows a 1- dimensional symmetry algebra for Z classes of interactions 
Fnm- The algebras and interaction functions can be represented as follows: 



: X = dt + 



(8) 
(9) 



nm ^nm^Unm i 



iUnmT'"' 



Fnm = Un,nfn,n{t, ^pq) , S.pq = , , (p, £ T \ {{n, m)}. 

(Upq) 

^1,3: X = Xnm{t) du„^ , 



J? — 



A. 



\q], {p,q) e r\ {(n,m)}. 



(10) 

(11) 

(12) 
(13) 
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Proof. We suppose that the system ([T]) has at least a one-dimensional symmetry group generated 
by vector field of the form Using the allowed transformations, we transform X into three 
inequivalent classes: 

a) For the case r(i) ^ 0, we choose the functions t{t) and Qnm{t) so as to transform T{t) 1 and 
\im{t) ^ 0. More precisely, we look for functions t{t) and Qnm{t) satisfying the following 
ODEs 

rit) t = 1, rQnm - KLm{i) ^ 2 ^ "iim^ Qnm{t) = 0. 

Hence, under allowed transformations, the vector field (O with T{t) 7^ is given by ([5]). We 
can now solve the remaining equation ([7]), for t = 1 and Xnm = 0, by applying the method 
of characteristics and we find function ([9]). 

b) For the case r = and a„„i ^ we choose the function Qnm(t) such that Xnm{t) — > 0, i.e. 



^nm^ nm \ 



.(t) + A„™(t) = 0. 

The vector field ([6]) is then given p^ . The remaining equation ([7]) gives (|lip . 

c) Finally, when t — and a„m = we already find the vector field (|12p and the remaining 
equation ([7]) gives us p^ . ■ 

Let us notice that vector field can be simplified in some cases. Indeed, the vector field 

can be transformed as X„m{t) A„„i(t) = p/^P,^^jA„„i(t) from the allowed transformations ([8]). 
Hence, if Xnm{t) is separable in terms of the discrete variables m and the continuous variable t, 
then we can transform Xnmit) into 1. 

We observe that the existence of a one-dimensional symmetry algebra restricts the interaction 
Fnm to arbitrary functions of 7 variables, rather than 8 variables in the original equation ([T]). In 
the next sections we will assume that the interaction Fnm, and one of the symmetry generators 
is already in "canonical form", i.e. they have the form ([T0| or (fT2|) with the corresponding 
interaction. We will illustrate how the interaction is further restricted by the existence of higher 
dimensional symmetry algebra. 

4 Abelian symmetry algebras 

Theorem 4.1. Equation ([7]) allows a 2-dimensional abelian symmetry algebra for 4 classes of 
interactions Fnm- The algebras and interaction functions can be represented as follows: 

M,\ ■ Xi ^ dt+ al^„\wnm9„„,„, X2 = anm^nmdu„^, 

Fnm = Unmfnmiipq), S.pq = (2) exp (l?[a(,^2,, a':^^\t) , {p, g) £ F \ {(n, m)}. 

{Unmr"" 

^2,2 : Xi ^ dt + anmUnmdunmi ^2 = , 
Fnm = am„Umn + e"'''"'* fnm{£,pq) , ^pq = P [u„„e"°""* , Ip,], (p, G F \ {{n, m)} . 
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Fnm — '^nrn fnm £,pq): {P: <?) ^ T \ {{n,m), [u + 1, m)}, 

Km ^ fnm{t,£,pq), £,pq ^ 'D[Un,nApq], (P, ?) ^ T \ { (tI, m) } . 

Two cases are not listed in these 2-dimensional symmetry algebras. One case corresponds 
to interaction Fnm with Xi and X2 of A4 4 (see above in the list of 4-dimensional symmetry 
algebras) and is not listed here since the symmetry algebra can be of dimension 4 with the same 
interaction, i.e. the algebra is not "maximal" . The second case corresponds to the degenerate case 
'^[^nm, ^n+i7n\ = ^ ^ith Xi and X2 of ^4 4. The generators are then given by Xi — Xnm{t)du„„, 
and X2 — Jm{t) Kim{t)du„,„ with 7,„ ^ 7m+i and 7,„. Hence, we obtain a non-isotropic system. 

Theorem 4.2. Equation (Qp allows a S- dimensional abelian symmetry algebra for 3 classes of 
interactions Fnm- The algebras and interaction functions can be represented as follows: 

As,! : Xi = dt+ ail^iUnmdu,,^, X2 ai^ijUnm X3 = anm'^nmdu^^, 
= Unm/nm {^pq ) , (P, (?) G T \ { (n, m) , (n + 1 , m) } , 
^pq = K„0"[''----k-n+l™)-"f'^--"lK,)"[''- exp (^V\o^^o!^l,,n,0^]^ ■ 

^3,2 : Xi= dt+ anmUnmdu^^, X2 = e""™*(9„,,^^ , X^ = Knme''"™'9„„^ , 
Fnm = aLi"nm + ^""""^ fnm{ipq) , {P, q) E T \ {(u, m), (n + 1, ITl)} 
^pq Ty\Unm^ j ^n+lm: ^pq\: ^nm 0; ^nm 7^ ^n+lmi ^nm 7^ 

^0,0 • — '-''nm nm^Unm ^ ^^2 — ^7im ^^m^Unm ^ — ^nm nm^Unmi 

Fnm = Unmfnm{t,^pq), (p, q) £ T \ {{n,m), {n + l,m), (n, m + 1)}, 

■T-^r (1) (2) (3) 1 

Again here, four 3-dimensional Lie algebras are not listed below. Three of them are not listed 
because they are not "maximal". The corresponding maximal algebras being given by Aq^^, ^4^4 
and ^4 5 above. The other algebra corresponds to the degenerate case Xn+im = A„m of ^4 5 
(without X4). The generators are then given by Xi = du„^, X2 = tdu„^ and X3 = Xm{t)du^^ 
with Am+i ^ A„i. Hence, we obtain a non-isotropic system. 

The discrete function Knm in ^3,2 depends on n and to, otherwise we obtain a decoupled 
non-isotropic system. 
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Theorem 4.3. Equation (QP allows a A- dimensional abelian symmetry algebra for 5 classes of 
interactions Fnm- The algebras and interaction functions can be represented as follows: 

^4,1 : Xi = dt + a,nmUnmdu„^, Xt = almUnm9„„^ , i = 2, 3, 4 

Km = Unrrifnmi^pq), (p, q) \ {{n,m), {n + l,m), {n,m+ 1)}, 

X M-^l^M.^^'^^U^ exp (-P[aW , 41+1, ■ 

^4,2 : Xi^dt + anmUnmdu^^, ^2 = e"""*9„„„, X,+2 = K^^e"'""*^,,^^ , i = 1, 2 

^lim = a^„„Um„ + e°-'"^*fn,n{^pq), {p,q) G r \ {(n,TO), (n + l,r7i), (n, m + 1)}, 

^4,3 • -^i = 0,nmUnmdu„,„ i i = 1, • • . , 4 

^nm = Un„ifnm{t,^pq), {p,q) G { (n - l,m), {n,m- 1), (n + l,m - 1)}, 

2 

^4,4 : X, - A(i(i)9„„„ , i = 1, 2; % ^ [Y. ""^i 5u,„„ , A: = 1, 2 

A^',)], (P, g) e r \ m), (n + l, m)}, I?[AW , Alfji J ^ 0. 

2 

wi/i ^ [^kjXlil + 2^kj\lil) - 0, det(cOfe,) ^ 0. 
^4,5 : Xi = au„„ , X2 = i9„„„, ^3 = A„m(t)9„„„, X4 = (cj(i)A„™(t) + (T(t))9„„^, 





1 






X(i) 




\(2) 


X(2) 


X(2) 
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1 / 1 /■* ct(s) \ 

\ ^ Jo V'^l'*) / 

A,i+i„i 7^ A„m, (p, e r \ {(n,™), (n + l,m)}. 

Theorem 4.4. Equation (QP allows a dimensional abelian symmetry algebra for 3 classes of 
interactions Fnm- The algebras and interaction functions can be represented as follows: 

A^.i ■ Xi ^ dt + ail^Unmdu^,^ , Xi = al'J„u«m9u„,„ , « = 2, . . . , 5 

Fnm = Unmfnrn{£.pq), {P, Q) G {("- - («, 771 - 1), (n + 1,771 - 1)}, 

)C[4^i,ai^ii„,ai''i^i,a^f!i„ + J f _pr„(l) „(2) (3) (4) (5), A 

A,2 : = a* + a„„u„„a„„„, ^2 = e"""*5„„„, ^ Kl*^e''"™*a„„„ , i = 1, 2, 3 

-Fnm = a^^U„™ + e"'"^"'*^fnm{Cpq), {P,^) ^ {{n ^ 1,"7), {n,m~ 1), (n + 1,771 - 1)}, 
C —T)\it p-i„„t (1) (2) (3) -, 1 

^5,3 : -'^i = alm"mna„„„., 1 = 1, . . . , 5 
-Film = Unynfnin{t,£,pq), {p,q) G {(71,771 - 1), (n + 1, 771 - 1)}, 

7-,r„(l) „(2) „(3) (4) (5), -nr„(l) „(3) „(<!) „(5)l 

(1) (2) (3) (4) (5), (1) (2) (3) (4) (5)1 



■^yanm+l) \an~lm+l) 

N-Pfa(i' Q<^' a'^' a<''> a'^'l / Nl3[a(i' a<^> a' 



Theorem 4.5. Equation (Qp allows a Q-dimensional abelian symmetry algebra for 5 classes of 
interactions F„m- The algebras and interaction functions can he represented as follows: 

Ae^i : Xi ^ dt+ al^^Unmdu^,^, Xi = a^^l^Unrndu^,^, 1 = 2, . . . , 6 

Fnm = Unmfnmi^pq), iP,q) € {(71,771 - 1), (71+1,711- 1)} 

/ \I7[a(2' a<^> a<''> a<^> aC^'l / N-I^fa'^' a<^> a<''> a<^> a(8)l 

X(u„„+i)'''°'""''''.+ i'"''''-i'"+i''''-i'"'''p5<(7t„_l„,+l) ^K„„.«„+i™,'i„™+i,'i„-i™:'^p,l 

/ ^■D^a<2) a<^' a<''' a<^> qC^'I/ ^-V^a'•^> a<-'> a'''> q<=' a'"' 1 



X exp 



y '^["'nm, «„+lm' "nm+1' "n-lm+1' "n-lm' "P9 I'' J ■ 



9 



^6,2 : Xi=dt+ anmUnmdunm . ^2 = e°""^*du^^ , Xi+2 = Ki*ie"'""*5„„„ , i = 1, . . . , 4 
Fnm = almUnm + e'''"^*fnm{^pq), ip,q) € {(n, m - 1), (n + 1, m - 1)}, 

<r _-nL, p-a™f (1) (2) (3) (4) , , 

^6,3 ■ Xi — ^nrn^nm^Unm ' ^ — ^7 • • ■ ) -^nm ^-^'nm./nm (^7 



ti+l''^Ti-lTTi + l'"-n-lm''^TiTTi-l'"'n + lm-lJ 



^ -Pfati' a'^' a'^' a'""' a'*^' a*''' 1 / ^ PFa'^' a'^' a'^' a<*> a<^' a"" 1 

sy(„, \ ■^l"nm><*n + lm''*n-lm+l''*n-lm'''nm-l'''n + lm-lJ \ ^l"nm ""n + lm '"nm + l ''^n- Im '"nm-1 '"n + lm-lJ 

'^\"'nm+l) {""n—lm+l) 

^(n. \ '^L'*nm'"n+lm'''nm+l'"n-lm+l'''nm-l''*n + lm-ll /„, \ l"nm '"n+lm '"nm+l '"n- lm+1 '"n- Im '"n + lm-lJ 



X (w„+1to-i) 



(1) (2) (3) (4) 



^6,4: Xi = Xl:l{t)du^^, i = 1,2,3; %= (j2^kdt)Xlil{t))du„^, k = 1,2,3 



1 + 





1 






X(2) 


j;(2) 


X(3) 


j;(3) 



1 ^n+lm 
x(l) 

71+ Im nm+1 

x(2) x(2) 

n+lm nm+1 

\(3) x(3) 



Cpq =T>[unm,>^n+imi^nm+i^^pq]^ ip^l) € T \ {(n, m), (n + 1, m), (n, m + 1)}, 

T^[Xnrn, ^n+\mi ^nm+l] 7^ 0' '^['^nm, '^I+Ito] 7^ 0, 

3 

mtt ^ (uJkjX^^l + 2d;fe,AW ) = 0, det{6jkj) + 0. 

2 



^n+lm ^nm \^n+lm ^nm ) ' ^{^nm, ^n+lm^ ^nm+l\ 
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i;(2) 

^nra 


X(2) 


X(2) 





1 


1 



' ~ T^[Unrm •^i-f-l,jn' '^im+li ^P?]' 



2 



with 

i=i 

Theorem 4.6. Equation (QP allows a 7-dimensional abelian symmetry algebra for 2 classes of 
interactions Enm- The algebras and interaction functions can be represented as follows: 

M.\ ■ Xi = dt+ al^„\wnm9„„,„ , Xi = a^^l^Unmdu„,^ , i = 2,...,7 

-^nm ^7unfnm{^^: 



'"nm-l'"n + l,m-ll 



/ N-I7Ia(^' a<^> a<''' a<^' a'*^' a'''' 1^ NP[a(^' a<^> o*"*' a<^> a'"' a<^' 1 

:(Unm+l) ^'''"'"'""+i"'''"-i"+i'°"-i"'°"™-i'°" + i-"«-i'(w„„i,„+i)^'''""'''"+i'"'''"™ + i'°"-i"'°"'"-i'''"+i."-i' 



^(Wn— Im) ^ 



a+lm''^7^m + l'"•7^-l7T^^-l 



(2) „W ^(5) ^(6) (7) , 



a'*"' 



X exD f-Pfa^'^ a*'^ a^^^ a^^' a*^^ a^^^ a'^^ u") 

A exp l^[Unrn, "„+lm' "nm+li "n-lm+1' "n-lm' "rim-li "n+l,™-!]''^ • 



^7,2 : = (9* + Onm-UnmSu^^ , X2 = e''""*a„„„ , = , i = 1, . . . , 5 

F -„2 , a„„.t^ -(i) _Q H) I (i) [i) , (i) 

^ _ „r -a„„i (1) (2) (3) (4) (5) , 1 

Theorem 4.7. Equation (Qp allows a ^-dimensional abelian symmetry algebra for 2 classes of 
interactions Enm- The algebras and interaction functions can be represented as follows: 



4 

^8,1 : X, = A,« (i)an_, i = 1, . . . ,4; = (^^fc, 

i=i 



fc = 1, . 





1 












n+lm 


a(i) 


A 


X(2) 


V(2) 
^nm 


x(2) 

n+lm 


x(2) 

nm+1 


A 


X(3) 


V(3) 
^nm 


n(3T 
n+lm 


X(3) 


A 


),(4) 


V(4) 


X(4) 


X(4) 
^nm+1 


A 



(1) 

n-] 
(2) 
n-] 
(3) 
n-] 
(4) 



nm , A„ 



(2) ;^(3) 



+ lm' nm+1' Vi— lm+1 



(4) 



+ fnni{t,£.pq), {p, q) E {{u - 1,to), (n, m - 1), {n+ 1,171- 1)}, 
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P -T)\n, X(2) ,(3) >(4), „r,(l) ,(2) ,(3) >(4) , ,„ 

2?[Aim, A|j|j^, A|f2i+l] 7^ 0, I?[Al2i, A|j|j,„] 7^ 0, 



with 



Aa,2 ■■ = 9„„,„, X2 = <a„„,„, X,+2 - AW^(i)a„_, z = 1,2,3 
3 

Ffe- (E^fc.WA(£W + f^^'Hi))<™, fc= 1,2,3 



A 



(1) 



A 



(1) 



A 



(1) 



I^2-2?KAi'ji„,A 



CX(1) _ A ,(3) -, 1 

n+lm '*n™ l^n+lm '^n™ 7 ' ■'^ L'^n™ i '^n+lm ; ^nm+l ' -"-"-Im+lJ 

2^3 • 2?[Aim, A|j|]^j„, l„m+l] • T^[u, ^nlim^ ^nm+lJ In-lm+l] 



\/^n+lm ^mn) ' -i^[^nm, A„^2mJ ' '-^V^rimT ^n+lmi ^nm+li ^n-lm+l\ 



,(1W2) 



l(l) n(2) 



,(3) 



I?3 • I?[u„m, An'ljlj^, l„m+lj 

TxW _ N r,r\(i) \(2) ^ ^ 

V^n+lm ■^nrn) ' -i^[Anm, ^n+lrm -"-nm+l 



= P[m„™, A|,^ji„, Aj,''^+i, Aj;'Zi„+i, Ipg], e {(n - l,m), (n, m - 1), (n + l,m- 1)}, 



>(2) 



,(3) 





X(i) 


















'D2 = 


X(2) 


X(2) 


\(2) 






X(3) 


X(3) 


X(3) 









1 


1 









1 


1 




n+lm / '^*nm' 


PfA^i) 


x(2) 1 1 


^0, 


2?[Ail^,A 


(2) 

n+lm ' 


A 1 


-Im+l] 7^ 



Theorem 4.8. Equation (1^ allows a 10 -dimensional abelian symmetry algebra for 2 classes of 
interactions Fnm- The algebras and interaction functions can be represented as follows: 



+0,1 : X, = \(:l{t)du^^, i = 



Yk = 



(14) 







aSS. 
aJS^ 


^,1 + 1 




Ai^i 


aJ."! 


4! 




A<fi 


aSI 





+1 




1 + 1 




1 






a'i> 




a"i) 






+1 


A„_lr 


n+l 




1 






A<2> 




a"=^) 




+1 


A„-lr 


n+l 




1 






A<3> 




a"^) 






+1 




n+l 




1 








a"") 




+1 




n+1 




1 




a(^> 




a"^) 






+1 




1+1 




1 



A^'^ A 



(3) 



,a: 



(4) 



n+lm ' nm+1 ' n — lm+1 ' n — Im 



,A, 



(5) 



(15) 
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= I'Km- Ai+lm,Ai^m+l.^i-lrn+l>^i-lm>49]. l) ^ {(n, TO - 1), (n + 1, m - 1)}, 

-r)r\(l) \(2) \(3) x(4) x(5) 1 / ^ x(2) >(3) ,(4) n 

'E'[AiTO, A^2t+l] 7^ 0' ©[Aram, A^^i^] ^ 0, 

5 

With J2 {^kjX^£ + 2d;fe,A« ) = 0, det{ukj) + 0. 
^10,2 : ^1 = i9„„„ , ^2 = ta„„„ , X,+2 = A£'^ (0^«r.™ > z = 1, . . . , 4 

4 

% = (E W) + '^'''(i))^., fc = 1, . . . , 4 



(16) 



^1) 



• 2?[Mnm, A 



(1) 



n+lm' lnm+1 



+ ©2- 



2^3 ■ 2^[Anm, A^_^]^„, l„m+l] 
2^[Anm) •^i+lm' ^nm+l] 



'^b^nmi -^K+lm' '^wm+l' In-lm+l] 

TxW _ ^ ^(3) ^ '] 

A^n+lrn ^nm; ' -i^lAnrnj ■^n+lm' J^n-lm+lJ 



^3 • P[A1to, A^j.^^, Inm+l] 2?2 ' P[Atoto, A^^^^, A^j^_)_^, 1 



(1) A (2) 



+ 2?4 



n — Im+lJ 



'^'[Alm, A^l]^^, Inm+l] 



+ 



T>\\<^^^ a(2) a(3) 1 , J 
■'^L'^nmj '^ri+lm' -"-n-lm+lj 

-n -nri(i) i^^) i i -nrx^^) x^^) ,(4) ^ 

• U[Anm-i '^n+lmi ^nm+l\ " -i^L^nm, ^n+lm' '^nm+l' ■'■n-lm+lj 



-Drx(l) X(2) 1 J -r)[\(l) ,(3) , 1 

■'^L'^nm) ^ra+lm' -"-nm+lj -t^L^nm) ^n+lm' ^n-lm+l\ 

T>U, X(l) X(2) >(3) -I 
-^["rim, ^n+lm' ^nm+li -^n-lm 



('x(^) _ x>-^M r)[x>-^'' x*-^-' x*-'''' x*-'*'' 1 

\^n+lm '^nmj J-^l'^nm, '^n+lm' ^n-lm+1' 



(1) N 7ir\(i) \(2) \(3) \(4) 





V(i) 




n+lm 




V(i) 








V(l) 




A^i) 

n+lm 


1^2 = 


j;(2) 


5,(2) 


X(2) 

^n+lm 


, 2^3 = 


V(3) 


X(3) 


X(3) 




V(4) 




^n+lm 







1 


1 







1 


1 







1 


1 



= 'DlUnm, >^nllm^ '^im+l ' ' '^n-lm' Ip?]' (P> ?) ^ {(n, TO - 1), (n + 1, TO - 1)}, 

^n+lm 7^ '^nm' ^[•^Im) '^i+lm' Inm+l] 7^ Oj ^[-^nm' '^i+lm' ^nrn+l' In-lm+l] 7^ 0, 

4 

'D[>^'nl, A^fji^, a22,+i, Ai'2i„+i, l„_i J ^ 0, with (d)fe,-A« +2d;,, A« j+cT^'^) = 0, det(d;,,) ^ 0. 



,(3) 



.(4) 
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Theorem 4.9. Equation (1^ allows a 12 -dimensional abelian symmetry algebra for 2 classes of 
interactions Fnm- The algebras and interaction functions can be represented as follows: 



Ai2,i : X, = A,W z = 1, . . . ,6; ^ (^^fcJ(^)A(f2^(^))5n_, k^l, 



,6 (17) 





1 










3;(i) 




^nm+l 


A 


\(2) 


\(2) 

^nra 


X(2) 
^n+lm 


X(2) 


A 


X(3) 


X(3) 


X(3) 


X(3) 


A 




3;(4) 


^n+lm 


^nm+1 


A 


n(5) 


V(5) 


X(5) 
^n+lm 


X(5) 


A 


\(6) 


X(6) 


\(6) 
n+lm 


\(6) 

nm+1 


A 



'^n— lm+1 

(1) 

n— lm+1 
(2) 

n— lm+1 
(3) 

n— lm+1 
(4) 

n— lm+1 
(5) 

n— lm+1 
(6) 

n— lm+1 



^n — Im 

n—lm 
\(2) 

n—lm 
\(3) 

n— Im 

n— Im 

A<^> 

n— Im 
\(6) 
n—lm 



''^nm— 1 

nm— 1 
\(2) 

nm— 1 
\(3) 

nm— 1 

nm— 1 
n(5) 

nm— 1 
\(6) 

nm— 1 



(18) 



fvw^'^^ A^^) a'^^ a^'^^ a'^^ a^^^ ]) \ f (t n 

f-nU, X(l) X(2) x(3) ^(4) ,(5) <(6) 1 

S ~ -^["n™' '^n+lm' ^nm+li '^n-lm+li ^n-lmi ^nm-li '^n+lm-lJi 



r)fX(l) x(3) ^(4) - (5) ^ (6) 1 , „ ^r, (1) ,(2) > (3) , (4) > (5) i , „ 

■^L^rimi ^n+lm' ^n-lm+1' ^n-lm^ ^nm-lJ 7" -^L^nmi ^n+lm' ^n-lmJ A 

r)fx(i) x(2) x(^) x*"') l=^n rjfxfi) x'^^ x^^^ l =^ n 7:)fx(i) x^^' l^n 

■'^[^nmi ^n+lmi ^n-lm+lJ 7^ ■'^L^nm' •^n+lm' ^rim+lJ 7^ ■'^L^nmi ^n+lmJ ^ 

6 

mt/i {u;k,\^± + 2c^fc,A« ) = 0, det(^fc,) ^ 0. (19) 



^12,2 : Xi = du„^, X2 — tdu„^, Xij^2 
5 



A^i(i)5n_, z = l,...,5 



(^a;fc,(^)A^2^(^)+'^fe(^))5n„„, fc=i,...,5 
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V(i) 



2?3 • V[Unrm ^i+lm' Inm+l] 



■nm — 

A, 



n+lm ^nm 



(^n+lm ^nm) 



('^i+lm -^nm) • 'Z^[AnTO, Aj^_ji;^^, 



(1) A (2) 



+ 



+ 



+ 



^3 • ^[Airrt, ^i+lnn Inm+l] 

r)r\(i) \(2) ^ ^ 
-i-'l.^nm) ^n+lm' -"-nm+lj 



77h/ A^l) 1 

-^["nm, ^n+lrriJ -^n 



-Im+lJ 



a^^) - X(l) ^ T)[X(1) X(2) ,(3) , n 

, V^n+lm ^ram; ' '-'[^nmi ^n+\m^ ^nm+l^ ^'"■-Irn+l] 



I'IAIto, A^i.j^^, l„m+i] 



(1) \(2) 



(4) 



7^r\(i) \(2) \W 1 
'^[^nmi ^n+lm' ^nm+\^ J-n-lm+l 



(3) 



2?3 • PfAiTO, A^^]^^, • 'D[\nm-, AJ^+Ito, A^^_|_i, 



(1) A (2) 



2?[Airn) A^^ijjj, Inm+l] " ^[Anmj A^_i!.i„, A),^_|_]^, \n-lm+l 



(1) a(2) 



v(3) 



_ • P[Anm, A^'^j.i^, Inm+l] _ I>2 ' P[A»V„ , >yn+lm- ^nin+\^ Iw-lm+l] \ 



^4) 



2-'[Arim) A^_^]^^, Inm+l] T^\^rni ^+lmi Anm+l' ln-lm+1 



(1) x(2) 



>(3) 



"Ds • V[\nm, A^+im; Inm+l] ' I'iAnm, ^\i+lm^ ^'nm+li ^n-lm+l] 



(1) x(2) 



■'^['^"mj ^n+lmi ^nm+l\ ' '^[■^nm, ^n+lnf -"-n-lm+l 



(1) \(2) 



v(3) 



l^[^nm, ^„+1to' -"-nm+lj 



1^5 



I?l • I?[Amn, A|f|j„, l„m+l] 



r)[x(i) x(2) 1 1 

i^[Anm, '^n+im,' ^nm+l\ 



Zl + fnm{t, 0> 



z = 



(1) ,(2) .(3) .(4) .(5) ^ 1 



T)U, x(l) x(2) x(3) 
^["nm; ^ra+lrrt' ^nw+l' ^n-lm+1' -'-n-lm 

_ ^ -nrx(i) x(2) \(3) 



lm+1' In-lmJ 
Inm— 1 1 



^1 = 



r)L, \(2) \(3) x(4) 

-^["wm; ^n+lrrn ^nm+li ^ra-lm+l' ^w-lm' 

-r)r\(i) \(2) \(3) \(4) 



77(1) r)rx(i) x'^) T^^) 7(4) 7(5) ^ ^' 

(3) 



Z2 — 2?[Ai"2ii Ai+i„, A^'2i+ii In-im+i] • ^'i, Z3 - I?[A^-2ii A*-^^ 



(1) Xl3j ^(5) , ,7 





V(l) 








V(l) 




A^i^ 

^n+lm 




V(l) 


A^i^ 

/»ri7n 




P2 = 


j;(2) 


X(2) 


\(2) 

^71+1771 


, 2?3 = 


\(3) 


\(3) 


x(3) 

^7l+l777 




\(4) 


X(4) 


X(4) 
^ra+lm 







1 


1 







1 


1 







1 


1 
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X(i) 








V(5) 


\(5) 


X(5) 







1 


1 



/ \(i) 7ir\(i) \(2) 1 1 ^ n 

T)[\(l) ,(3) , wn 7^r\(l) \(2) \(3) x(4) -, i / r, 

'^[■^nrm '^n+lrm ^nrn+lJ '-n-lm+l\ 7^ "j ■'^ [^nm 1 ^n+lm i ^nm+1 ' '^n- Irn+l ImJ 7^ "i 

5 

5 Nonsolvable symmetry algebras 

The vector fields of equation ^ should have the form ^ . We can now ask whether it is possible 
to obtain simple symmetry algebras from these vector fields. We obtain the following theorem. 

Theorem 5.1. Equation ([7]) allows only one simple Lie algebra, sl(2,M), given by: 
NSs^i : Xi = dt, X2 = tdt + ^u«m9„„„, X3 = t^dt + tunmdu^^, 



(20) 



F, 



' fnm{£,pq) 1 £,pq 



nm 3 Jnm\^pqji sp<3 



{P,q) G r\ {(n,m)}, 



[Xi,X2] — Xi, [XijXj] — 2X2, [X2,X3] — x^. 

We can now look for additional symmetries by considering nonsolvable symmetry algebras for 
equation ([1]). A nonsovable Lie algebra must contain a simple subalgebra, i.e. the Lie algebra 
sl(2,R) of 5*3,1 in our case. Therefore, we add new vector fields Yi of the form ([6]) to NSs^i. 
These vector fields {Yi} forming the radical of the new Lie algebras. The following theorems give 
the possible nonsovable symmetry Lie algebras for equation We only list the radical of the 
nonsolvable Lie algebras since all of these algebras have sl(2, M) of the form ((20|) as subalgebra. 

Theorem 5.2. Equation ([7|) allows a 4- dimensional nonsolvable symmetry algebra for 1 classe of 
interaction: 



nm — ^nm 



I3la„„,.l„+i„| 



fnm (£,pq ) j 



/ \a„+im t \-an 

^pq = (^/„+i™)^["— , (p, g) e F \ {{n, m),{n+ 1, m)}. 

Theorem 5.3. Equation ([7]) allows a 5 -dimensional nonsolvable symmetry algebra for 2 classes 
of interactions: 



F„, 



(1) 



(1) (2) — ; : 



y^nraj y^n+lrnj y^nm+1 j 

x/„m(Cpg), e r\ {(n,m), (n + l,m), (nm + 1)}, 

spq — \anm} \^n+lni) \'^nra+l) \^pq J 



NS5^2 ■ Yl — du„^, Y2 — tdu„^, Fnm — (T^iUrun, In+lm]) fnmiC, 



ml 



^[^nm; In+lm] 



, (P,?) e r\ {(n,m), (n+ l,m)}. 
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Theorem 5.4. Equation ([I]) allows a 6 -dimensional nonsolvable symmetry algebra for 1 classe of 
interaction: 



Fn.m. — Un 



nm , 



Tir (1) (2) (3) , ~r (1) (2) (3) 1 Tir (l) (2) (3) 

\an+lmj \anm+lj 



X (u„_lm+l) 



fnmi^pq) 7 



4: 

(1) (5) (5) {P,q) G {("-- 1,"^), {n,m- 1), (n + 1,771- 1)}. 



Theorem 5.5. Equation ([7]) allows a 7 -dimensional nonsolvable symmetry algebra for 2 classes 
of interactions: 

NSr.i ■.Yi = alll^Unmdu„r„, 7 = 1,..., 4; (p, g) G {(?7, tti - 1), (?7 + 1, 777 - 1)}, 



^ ^■D^a<l' a<^' a'^"' a'"' 1 1^ \-I?[a(i' a<^' a<^' a<*' 1 , 



4 



"™ r,r (1) (2) (3) (4) 

J-^[a-nm, O-n+lrm '^ri-lm+1' J-n-lmJ 

NSt,2 : = du^^, Y2 = tdu^^, — Knmdu„rnJ ^4 = Hnmtdu„ 



^ nm — y^i^nniT ^n+lrn: ^nm-\-l\) Jnmy^pq)^ Sj 



pq 



^[^nrrn ^n+lmi l^im+l] 
VlUnra, Kn+lni, Inm+l], (p, q) ^ T \ {(?7, 77l), (77 + 1, 777), (77, 777 + 1)}. 
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Theorem 5.6. Equation ([I]) allows a 8 -dimensional nonsolvable symmetry algebra for 1 classe of 
interaction: 

iV 5*8,1 : Yi = a^nliUn,ndu„m. I « = 1, • ■ • , 5; 



I \V\oS'-'> a<^> a<^' a**' a<^' 1/ ^-X)^a<l) a<^' a<^> a'"' a'''' 1 

/ ^X'^a<l> a<^' a<^' a<*' a'""' W \I7[a(i' o!^'^'^ o<^' a'"' a'^' 1 

X (u„_i„j)''^L""™'''"+i™'°"'"+i'''"-i'^+i'''"'"-i'(M„m_l)^'°"'"'''"+i™'''"™ + i'°"-i'"+i'°"-i'"' 

^fnm{S,), Ctnm = (i) [31 [51 (I) ~) ; 

T-,r (1) (2) (3) (4) (5) . , „r (1) (2) (3) (4) (5) , 

/ N-I5Ia(i) a<^> a<^> a<''> o<^' 1 i, ,1^ ^VU^^ a<^> a<^> a*"' a'''' 1 a., 



(5) 

xK+i,„-i)-''"'"-'''"+i"''^-+i''^"-i'"+i''^"-i'' 



Theorem 5.7. Equation ([7]) allows a 9 -dimensional nonsolvable symmetry algebra for 1 classe of 
interaction: 

NSg^l : Yi = du,^„^, Y2 = tdu„^, I3 = K-nmdu,^„^, Y4 = Knmtdu„^, 

F — 



f'nl (1) (2) -, A f — -^Prtm, ^n+lm; ^? 

I i^[U„,„, i„_l„i-|_lj I Jnm[C,pq), ?P9 — (j^) (2) 

2?[u„m, '^n+lm; '^nm+li ln-1 



+i> ^ pgJ 

Im+lJ 

e {(n- l,m),(n,rn- l),{n + l,m- 1)}. 



Theorem 5.8. Equation ([7]) allows a 11 -dimensional nonsolvable symmetry algebra for 1 classe 
of interaction: 

NSii^l : Yi = 9ti„„, Y2 — tdu„^, Y^ = Knmdu„^, Y4 — nimtdu^^, Y^ = Knmdu„^, 

Yq — l^nmtdunrri^ ^^m^Unm^ ^^mtdu^^^ 

(21) 

^nm ~ ^^ *=li2,3; Fnm ~ [^nm J '*n+lm ' '*nm+l ' '^n- Im+l ' 1"-!™] ^ fnm{£,pq\ (22) 

7-)r (1) (2) (3) -, 1 

E - ^^^"'"'^"+lm.^nm+l:^»-lm+l:-^pd f» C / fn m - 1 Wn + 1 m - 1)1 

spg — (2) (3) i' \Pt h) ^ IVh "I- i^), \n -t -L, III i)}. 

'-^[Unrm '^nm+l' '^ri-lm+1' -'^n-lmj 
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Theorem 5.9. Equation ([7]) allows a 13 -dimensional nonsolvable symmetry algebra for 1 classe 
of interaction: 



Yq — t^^rntdunm^ ^ ~ ^^nrn^^Unrn^ Y^ — t^^mtdu^^^ Yf) — i^nrn^UnTni Y\{) — i^rimt()\ 



(23) 



r)r (1) (2) (3) (4) , , 



K^*) =0 7 = 1 4- f = 

"'nm ^^ l — ±, . . . , 'l, <; — 



V\v K^^) K^^^ K^''^ 1 J 



6 Conclusions 

Group theoretical methods have been used to classify equation ([T]) according to their symmetry 
groups. Abelian and nonsolvable symmetry algebras C have been considered; the class of linear 
equations ([T]) has been excluded. The results of the symmetry classification are summed up in the 
following table: 



dim>C 


Abelian 


Nonsolvable 


Total 


1 


3 





3 


2 


4 





4 


3 


3 


1 


4 


4 


5 


1 


6 


5 


3 


2 


5 


6 


5 


1 


6 


7 


2 


2 


4 


8 


2 


1 


3 


9 





1 


1 


10 


2 





2 


11 





1 


1 


12 


2 





2 


13 





1 


1 



Table 1. Results of the symmetry classification of equation ([T]). 

Classification given in this paper has applications in solid-state physics. The interest in two- 
dimensional systems has its origin in possible applications to magnetic systems as well as to ab- 
sorbed layers. The theoretical models treated so far are mainly spin-models and range from Ising 
systems with competing interactions to planar Heisenberg-models [501 HI] ■ Models involving equa- 
tion of the form ((!]) appear in references [6l [71 [8] where analytic and numeric calculations were 
performed. Lie point symmetries obtained in this work could be considered to obtain analytical 
solutions of these models by using symmetries to generate new solutions from a known one or by 
using the 'symmetry reduction method'. Moreover, some interactions found in this work could be 
considered as models with appropriate symmetries. 

A continuation of this study is in progress. We know that the existence of many symmetries is 
an indication of integrability. Consequently we can ask ourselves which of the equations that are 



19 



completely specified by their Lie algebras, and therefore that have many symmetries, are integrable. 
Completely specified equations to be considered for abelian and nonsolvable Lie algebras are, 
respectively, the following: 

^6,3, ^7,1: ^7,2, ^12,1 and Ai2,2, 

and 

NSs,!, and NSi3,i. 

Finally, a further task is to complete the classification, that is to treat the nilpotent and the 
solvable Lie algebras. 

Appendix 

This appendix includes the details of the proof for one abelian algebra and one nonsolvable algebra 
in the higher dimensional cases. 

For the abelian Lie algebra we are considering the proof of Ai2,i of theorem 14.91 Since the 
procedure to obtain this classification is to proceed by dimension (for each type of algebras: abelian 
or nonsolvable), we suppose that we already have obtained the algebra ^io,i of theorem 14.81 We 
add a new vector field of the form i.e. Z = T{t)dt + [(§ + a„„i) + A„m(t)] to the 

symmetry algebra (fH)) . 

Considering the commutation relations [X^, Z] — and [Y^, Z] = for i = 1, . . . , 5 we obtain 

5 5 

tAII = A^'^(^+a„,„) and r ( ^ A^^+Wfcj aI^],) = ( ^ A^^^i) +a"m) • (25) 

We separate the proof in two cases: 

A) The case t = 0. From the first preceding equations we easily find Z = \'i%m{t)du„,^, where 
Anm ■— Aim(t) for convcnicncc. We want now to solve the remaining determining equation 

(p,9)er 

where Fnm is the interaction psp of ^lo.i. This equation is equivalent to 

r,r\(l) 1,(2) x(3) i(4) ,(5) 1 

^ ^['^nmi '^i+lm' '^i-lm+li '^ri-lmi '^pg'']^Cpg/(^J Cpg) (26) 

(p,g)er' 

-X- T)rx(l) X(2) x(3) \(4) A5) ,(6) i _ „ 

where F' := {(n, m — 1), (n + 1, m — 1)}. Using then the method of characteristics we find 
interaction (jlSp with symmetry algebra of dimension 11: where the vector field Z := ~ 
^nmdu„^ is added to the Lie algebra of ^io,i- 

B) The case r 7^ 0. By the allowed transformations we choose t such that rt — 1 which implies 
T = 1. Moreover, we choose Qnm{t) for which anmQnm.it) + ^nm{t) — Qnm = and we obtain 
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Z — dt+anmUnmdu„rn ■ FiTst equation of ((25|) implies that Xnm — Knrn.e"'"'"^ i = 1, . . . ,b where 
Knm is an arbitrary function of n, m. Again here, using Pnm in the allowed transformations, 
we can normalize Knm to 1. Therefore, the Lie algebra we are considering contains the 
subalgebra Aq^2 and the classification has been already done in this lower dimension. 

We are now looking if, for the interaction (fT5)) . an additional vector field Z = T{t)dt + [(§ + 
+ Anm(i)]9u„„ can be added to the symmetry algebras obtained in the case A. The 
calculations are similar to those presented above. The commutation relations considered are then 
[Xi. Z] ~ Q and [Yj , Z] = for i = 1, . . . , 6 and j — 1, ... ,5. Remaining equation ^ implies: 

7-,rx(l) x(2) x(3) ,(4) ,(5) ,(6) n 

■t^[Arijn, ^n+lnn ^nm+li ^nm-l\ 

^ '^[■^nmi ^n+lm' ■^nm+l' \i-lm+l' ■^n-lm' ■^nm-l^ ^n+lm-li'^iJ y'-' 'i) ' J 

{p,q)er' 

-X- T)[\(l) ^(3) ,(4) , (5) ,(6) ,(7) n _ „ 

where Xnm '■= Xnm{t) and Fnm is given by interaction (|18p . Interaction (|18|) is invariant when Z 



is added to algebra of A if Xnmit) = i^6jit)^nm{t). Equation fI7^ then becomes equivalent 



to condition p9|) of ^12,1. Hence, in this case the vector field Z :— Yq = {J2j=i ^6j^nrn)du„^ is 
added to the vector field of the case A with interaction (fT5)) . □ 

We are now considering the details of the proof for the highest dimensional nonsolvable case, 
i.e. for theorem 15.91 We suppose that we already know classification for dimension 11, i.e. 
theorem 15.81 has already been shown. We add a new vector field of the form Yg — T{t)dt + 
+ Xnm{t)] du„^ to the vector fields ([21]). The Levi theorem [lH [23] tells us that 
every finite-dimensional Lie algebra £ is a semidirect sum of a semisimple Lie algebra S and a 
solvable ideal (the radical R): 

£ = 5 C> i?, [S,S]= S, [S, R] C i?, [R, R] c R, 

such that we have 
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[X,,Y^] = Y,a^kYk, i = 1,2,3 and [f„ fg] = ^ /S^feFfe, j = l,...,8 
fc=i fe=i 

where aik,Pjk are real constants. The commutation relation [^1,19] gives us that 



r(t) = Toe"i»*, aiQttnm = 0, A„„(<) 
where 



tile^is* - 4- (aigA^nl + aiQB^nlt + B^nl) , ai9 ^ 0, 

lR(l)f2, 4(1)^ , ^(4) „ _n 

2^nml' "T ^nm'- ~r Knm, "ig — U, 



/((») _ , (1) , (2) , (3) 

Anm '■— Oiii + Oii^Knm + Ui^Knm + CXilKnm, 

Tj{i) , (1) , (2) , (3) 

£)„m :— + CXi^Knm + Q!i6'*nm + CXisKnm, 

for i = 1,2,3 (functions with i — 2,3 will also appear in what follows) and Knt™ is an arbitrary 
function of n, m. Considering now the commutation relation [X2, Yg] we obtain 

aigTo = 0, (a29 + 1)to = 0, a29anm = 0, 
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and 



K^^^ -0 

(2a29 - l)Bi\l - 2aigBl?i = for aig ^ 0, 

(1 + 2a29)ai9A^'^ + (1 + 2a29)Bi!i + 2algA';^l = 

(|-a29)(Bi'2)2 = 

("29+i)«""i + ^nm for aig = 0. 

(a29 - ^)Al^l + B^i = 0. 
From the commutation relatfon [X3, 19] we find 

To = 0, asQanm = 



and 



(1 + ai9a39)s4U - afgSi^js + aig^l^^ = 



for q:i9 0, 



B^'^ - 



nm 

Knm + Sim + ^39^1™ = for Q!19 = 0. 

Anm + 0:39 Knm = 

Finally, for j = 1, . . . , 8 the commutation relations K , ig] imply pjganm = and 
I £>^^'^ +d oA'^'^ -0 tor J -1,3, 5, 7 



for j = 2, 4, 6, 8 



where 



nW) ._/?,/? ^(1) I ^(2) , ,3 ^(3) 



(0) 

nm • — 

Since A^m, -Bnm, C'^m atiA^Dnm are linear combinations of linearly independent functions 
appearing in the vector fields Yj, we can use linear combinations to simplify Yg. In the generic 
case, i.e. for any aig, we have that a„„i = Cnm and using linear combinations wc can transform 
dnm to zero. Moreover, we have that r = such that Fg = ^nm{t)dur^mi where Xnm depends 
on aig. In the case aig ^ we can transform Anm = BnX = by linear combinations. From 
equations obtained previously, this implies that Anm ~ Bnm = for i = 1, 2, 3. Since K^nm ~ 0, no 
additional symmetry is possible when aig 7^ 0. When aig = 0, one equation gives us that Bnm = 
and by linear combinations we can transform Anm to zero. Therefore, we find 



and Knm ■= 1 for convenience. 
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Note that we cannot use allowed transformations to simplify Yg. All allowed transformations have 
n already used to simplify vector fields, in 
From the remaining equation ([7]) we have 



been already used to simplify vector fields, in particular to obtain Knm — 1 in Yi and Y2 



- V K^'^^d F 
(p,'z)er 

where Fnm is given by (|22p . Using the method of characteristics we find that the new invariant 
function when Yg is added to iV 5*11,1 is given by ((24)l . 

Let us now verify that the Lie algebra {Xi, X2, X3, Yi, . . . , I9} is 'maximal' or not, i.e. if we 
can add another vector field of the form Yio = T{t)dt + [(§ + anm) + ^nm{t)\ with the 

same invariant function (1241). From the Levi theorem, we have 
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[^.,?io] = ^5,fe?fe, i = l,2,3 and [Yj.Y^] = Y.I^A 3^^,- 

k=l k=l 



where aik,Pjk are real constants. Now defining 



7(^) 


:= an - 












■= - 


L - (1) 


I ~ (2) 
r CXiQKnm 








■■= 








1- ^("'^ 
r Pjgn-nm 















the commutation relations [X^, l^io] and [Yj , Yio] give us the same set of equations obtain previously 
replacing A^iln > ylnm, . . . , Dnm *■ DnmjOiik > osik, Pjk * Pjk and Knm * '^rim for i = 1,2,3, 
j = 1, . . . , 9 and k = 1, . . . , 10. Again here we have that t = and a„m = Cnm- Using linear 
combinations we can transform to zero such that Yio = ^nm(t)dunm- For ai lo = we can 
transform A„m to zero by linear combination. In the case ai lo 7^ we can transform Anm to zero 
by linear combinations but not Bnm (since this function does not depend on Kmn)- Therefore, we 
obtain 

Yw = n'nltdu^^ ■ 

The remaining equation ^ 

with Fnm given by is identically zero. This complete the proof for the nonsolvable case of 

theorem [CTl □ 
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